UTPT-97-12 

DOE/ER/41014-21-N97 

hep-ph/9707313 

Power Counting in Dimensionally Regularized NRQCD 



Michael Luke" and Martin J. Savage'' 

(a) Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7 
(h) Department of Physics, University of Washington, Seattle, WA 98195-1560 U.S.A. 

(July 1997) 



Abstract 

We present a scheme for calculating in NRQCD with consistent power 
counting in the heavy quark velocity v. As an example, we perform the 
systematic matching of an external current onto NRQCD at subleading 
order in f , a calculation relevant for the process e^e~ —>■ hadrons near 
threshold. Consistent velocity power counting in dimensional regulariza- 
tion is achieved by including two distinct gluon fields, one corresponding 
to gluon radiation and one corresponding to an instantaneous potential. 
In this scheme power counting is manifest in any gauge, and also holds 
for non-gauge interactions. The matching conditions for an external vec- 
tor current in NRQCD are calculated to 0{g'^v'^) and the cancellation 
of infrared divergences in the matching conditions is shown to require 
both gluon fields. Some subtleties arising in the matching conditions at 
subleading order are addressed. 
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I. INTRODUCTION 



Non-relativistic QCD (NRQCD) ||T|] is a powerful tool for analyzing the dynamics of 
systems with two or more heavy quarks at momentum transfers much less than their mass. 
Such systems are more complicated than single-heavy quark systems described by the heavy 
quark effective theory (HQET) because the quarks scatter via the QCD potential. This 
introduces infrared divergences in heavy quark scattering near threshold in HQET which 
must be regulated by resumming an infinite number of insertions of the kinetic energy 
operator. Since this operator is subleading in l/mg, this violates HQET power counting. 
This kinematic regime of QCD is of interest for a number of physical systems, including 
quarkonium, e^e~ — > hadrons near threshold and nonrelativistic QCD sum rules Similar 
techniques are also of interest for other nonrelativistic systems, such as positronium [^] and 
low-energy nucleon-nucleon scattering |p. 

A concrete example where HQET power counting fails is provided by an external current 
in QCD coupled to gFg, where F is some Dirac matrix. In processes with a single incoming 
and outgoing heavy quark (such as semileptonic b ^ c decay) there is no potential scattering 
and HQET is the appropriate low-energy theory. Loop graphs such as Fig. |l](a) are well- 
defined in HQET, and the matching conditions for this current in HQET are currently 
known to 0{as, l/niQ) 0. On the other hand, the one-loop correction to quark- ant iquark 
production by the same current near threshold cannot be correctly described in HQET; the 
one-loop graph in Fig. ^(b) is infinite when the four-velocities of the quarks are the same, 
and gives rise to the well-known infinite complex anomalous dimension for the current [0. 
The appropriate low-energy theory for the second process is NRQCD, which treats potential 





(a) (b) 
FIG. 1. Two of the one loop contributions to the matching of an external current in HQET. 

scattering near threshold properly. 

In this paper we consider the general problem of matching an external fermion- 
antifermion production current in a nonrelativistic theory. We pay particular attention 
to power counting in the NRQCD expansion parameter v, the relative three- velocity of the 
heavy quarks. Power counting in NRQCD is less transparent than in HQET, for several 
reasons. First, since v is not a parameter in the Lagrangian, power counting is not manifest 
in the NRQCD Lagrangian, although there are simple rules for determining the v scaling of 
an operator f^. Second, the power counting for on-shell gluons differs from that of virtual 
gluons contributing to potential scattering pHlII and in order to have simple v counting 



this distinction must be implemented at the level of the Lagrangian. Also, in order to retain 
simple V counting beyond tree level the theory must be regulated with a mass-independent 
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regulator such as dimensional regularization, instead of the usual momentum cutoff (other- 
wise divergent loop integrals change the power counting of Feynman graphs in the effective 
theory). 

The paper is organized as follows. In Section |I| we discuss velocity power counting in 
NRQCD, and the relation between the results of Refs. and In order to maintain 
manifest v power counting in dimensionally regulated NRQCD we introduce two sets of gluon 
fields, corresponding to propagating gluons and an instantaneous potential. In Section ^ 
we discuss the matching of an external current in nonrelativistic Yukawa theory (NRY), and 
show that the dependence on the infrared regulator vanishes in the matching when both 
gluon fields are included in the low-energy theory. We consider this theory both because it 
is simpler than QCD, as well as to stress that manifest velocity power counting does not 
depend on working in any particular gauge, such as Coulomb gauge. In Section ^ we match 
an external current in NRQCD to 0{g'^v'^), and explicitly show that the low-energy theory 
reproduces the nonanalytic behaviour of QCD to this order. Finally, in Section |V] we present 
our conclusions. 



II. VELOCITY POWER COUNTING 

In NRQCD, the power counting of terms in the Lagrangian is different from HQET, al- 
lowing potential scattering near threshold to be correctly described. Operators are classified 
according to how they scale with the three- velocity v instead of l/rriQ 0. Since the kinetic 
energy of a nonrelativistic state is proportional to v"^ while the momentum is proportional 
to V, space and time derivatives scale differently with v, and power counting is not manifest 
in the NRQCD Lagrangian, 

jC- = 4 (^zDo + ^D^^ tfjf, - -G^'^^G^.^a + -Cg.f. + . . . . (2.1) 

Nonrelativistic fields are distinguished here from fields in the full theory by the subscript h. 
£g.f. is the gauge fixing term, Z)^ = is the gauge-covariant derivative {A^ = A^T'^) 

and the dots denote higher dimension operators whose matrix elements are suppressed by 
powers of v. 

In Ref. a rescaling of the fields and coordinates was introduced to make the v counting 
of operators manifest at the level of the Lagrangian. The NRQCD Lagrangian was written 
in terms of new coordinates X and T and new fields and A, where 

X = A,X, t = XtT, = Xq^h, A' = Aao^°, A' = XaA' , (2.2) 

and where A^; = l/rngv, \t = l/rnqv^, Xq = A~^/^ and A^ = A^o = {rngXl) In this 
form the v scaling of operators is manifest in the Lagrangian. In Lorentz gauge, 

-\ (d,A] - d,A\ - gV^fakcAUf + \ {^^Al - vd.At - gV^fabcA\Al 
-^{vd'A'', + &A^^' 
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-^[d^A^)\0{v,gV^). (2.3) 

Of course, there is no physics in a simple rescahng. However, in order to have an 
effective theory in which the v power counting is manifest, the additional prescription that 
terms which are subleading in v be treated as operator insertions must be added (otherwise, 
loop graphs evaluated in dimensional regularization will mix powers of v). As was noted 
in [jlO|, once this prescription is added, the rescaling in Eq. ( p.3[ ) misses important physics. 



While it provides the correct description of virtual gluon exchange corresponding to an 
instantaneous potential, it fails to correctly describe on-shell gluons. The problem is that 
the pole in the full gluon propagator occurs at /cq = k^, or, in terms of rescaled variables, 
v'^Kq = K^. Unless the time derivative in the gluon kinetic term is treated exactly instead 
of as an insertion (which would violate manifest v counting), amplitudes in the effective 
theory do not have the correct branch cut corresponding to physical gluon propagation, so 
the effective theory cannot describe on-shell gluons. 

In Ref. |Tl[] it was demonstrated that this problem is avoided by a further rescaling (in 



Coulomb gauge) of the space coordinates of only the transverse components of the gauge 
fields. In the language of Ref. fT^lQ) this corresponds to the rescaling 

X = XtY, t = XtT, A' = Xt^A\ (2.4) 

The kinetic term for the ^'s is canonically normalized, 

- J d'Y dT [^(9,^1 - djA^ - gfatcAlA'^y + ^ [doA]) 



-'kin 



2.5) 



while the transverse gluon-quark interaction may be expanded in terms of multipolesQ 



An* = -^(7t;^I(X,T)V,^,(X,T)A(^X,T) 

= -'-gv^iiX, T) V,vl/,(X, T)[l + vX-V + ...] MO, T). (2.6) 

Note that three-momentum is not conserved at the multipole interaction vertex, since the 
theory breaks translational invariance, although energy is conserved. Once the multipole 
expansion is performed, loop integrals in dimensional regularization do not change the v 
scaling of a graph determined by the vertices. 



^The Lagrangian in was written as an expansion in powers of 1/c instead of f ; however, the 
two descriptions are equivalent. Note that as c — > oo, = /A-kc — > 0, whereas as — > 0, 
remains constant, so factors of Ug scale differently with 1/c than with v. 

^At the level of Feynman diagrams, the observation that on-shell gauge fields couple via the 
multipole expansion was made in Ref. [^. 
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There is, however, a subtlety arising from the multipole expansion. Since three- 
momentum is not conserved at the vertices, transverse gluons cannot alter the three- 
momenta of the heavy quarks. Potential scattering via transverse gluon exchange therefore 
does not occur in the low energy theory. Since the amplitude for potential scattering is not 
analytic in the external momenta, it cannot be reproduced in NRQCD by the addition of 
local operators. Both potential scattering and real gluon emission are long-distance effects, 
so in order to correctly describe the infrared physics of QCD in the nonrelativistic limit the 
instantaneous potential due to transverse gluon exchange must be added to the effective 
theory. This may either be done by explicitly including spatially non-local operators in 
the nonrelativistic theory (the need for which was discussed in [|1T|), or by reintroducing a 
second gluon field which couples according to Eq. (|2^ ). 



The need for two distinct gluon fields is understood by comparing the energy and mo- 
mentum of radiated gluons compared to those involved in potential scattering. In the nonrel- 
ativistic regime of QCD there are nonanalytic contributions to scattering amplitudes arising 
from gluons in two separate kinematic regions, both with energy of order mv'^. The gluons 
with spatial momenta ~ mv are far off shell and contribute to potential scattering. In con- 
trast, the gluons with spatial momenta ~ mv'^ may be on-shell, describing real radiation, 
but do not contribute to the scattering of quarks with three-momenta of order mv (in the 
limit f — >■ 0). Each of the rescalings discussed above treats one of these kinematic regimes 
correctly, while missing the physics of the other. As argued in Ref. fl^, it is not possible 
to describe both regimes via any single rescaling. In order to correctly describe the infrared 
physics of QCD in the nonrelativistic limit two separate gluon fields must be included, which 
will be referred to in this work as "potential" and "radiation" gluons, corresponding to the 
two different kinematic regimes described above. 

In this approach the heavy quark Lagrangian in NRQCD in Lorentz gauge is therefore, 
in standard (unrescaled) units, 

^ = 4 (^^0 + ^ - gA% - gAliO, t)^ (v^A^p - A^p)' + \ {yA%f 

--G'rG,.r - ^ (V ■ Ap)^ - ^ {d.A^^f + 0{v, g^v) , (2.7) 

where the subscripts P and R denote potential and radiation gluons, G^^ is the field strength 
tensor for radiation gluons, and a is the gauge parameter. For practical calculations, this 
version of the Lagrangian is much more convenient than the rescaled version. The rescaled 
theory simply guarantees that loop graphs computed in the unrescaled theory, with the 
appropriate terms treated as operator insertions, will have v scaling determined by the 
vertices. 



III. YUKAWA THEORY 

Coulomb gauge is usually used in NRQCD because in this gauge Aq exchange corresponds 
to an instantaneous potential. Once v counting is performed as in the previous section, 
potential gluon exchange is instantaneous in all gauges for both transverse and longitudinal 
gluons. More generally, potential scattering proceeds via an instantaneous interaction even 
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in theories with no gauge freedom. To illustrate this in a theory which is simpler than QCD, 
in this section we consider a nonrelativistic Yukawa theory (NRY) of a massive fermi field 
coupled to a massless scalar 

As a warmup for our calculation of the matching conditions for e^e^ — > hadrons in 
QCD, we consider here the matching conditions for an external current coupling to tp'j'^ijj 
in NRY. In processes with a single incoming and outgoing fermion there is no potential 
scattering and the analog of HQET for a Yukawa interaction is the appropriate effective 
theory. However, as has been discussed, the 1/v behaviour of potential scattering near 
threshold cannot be correctly described in HQET. In NRY, the 1/v behaviour is reproduced 
by potential scalar exchange. There are also infrared divergences in the full theory due to soft 
scalar bremsstrahlung, which are reproduced in NRY by radiation scalars. In this section this 
is demonstrated explicitly. The theory is regulated in the ultraviolet by working in d = A — e 
dimensions, and in the infrared by introducing a small scalar mass m<^. (The theory could 
be regulated in the infrared with dimensional regularization, as at the end of this section, 
but this obscures the distinction between the infrared and ultraviolet divergences.) 

The Lagrangian in the full theory is 

C = ip{ip — m)ip + ^d'^ipdfj,ip — gipip^ , (3.1) 
while the nonrelativistic Yukawa theory (NRY), to leading order in the three- velocity v, is 



/ 1 1 

-^NR = i^do + —j + -d^^Rd^^^R - -(V^p)^ - gii)li)h^p - ^2V'iV'h¥'R(0, t) , (3.2) 



where <^p and (p^ are the potential and radiation scalars, respectively. There is a similar 
kinetic term for the anti-fermion field Xh- (Note that iph annihilates incoming particles, 
while Xh creates outgoing antiparticles) . At tree level gi = g2 = g- 
The external vector current in the full theory, 

J^tfj^^tfj , (3.3) 

matches onto a number of terms in the low energy theory,^ 

J, {xWA + ^Wxh) + Jo + xixh) + 0{g\ v""). (3.4) 

The one-loop matrix element of the current in the full theory is given by the diagrams 
in Fig. ^. Performing the loop integration, the infrared divergent part of the vertex graph is 

'-^'^ = '^^'"^'^^ ^(P^hMP2) J (t).(p_^2)(p + 2pV^)(F-2p,.A:) + • • • ' ^'-'^ 

where terms finite as have been neglected. Combining denominators in the usual 

way and performing the loop and Feynman parameter integrals, the infrared divergent term 
is found to be 



^Spinors in the full theory are normalized such that Ur{p)us{p) = —Vr{p)vs{p) = Srs- The 
two-component spinors are Uh^ = = (o)' '^h2 = ^fti ~ (i)- 
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FIG. 2. One loop diagrams in Yukawa theory. 



■.av_ r 



iA 



47^2 



Ji M7*f 



tanh ^ f ^ j lii^^-p + • • • 



where 



/3 



(3.6) 



(3.7) 



is the magnitude of the three-velocity of each fermion, and s = {pi + is the invariant 
mass of the fermion- ant if ermion pair. The infrared divergent piece of the wavefunction 
renormahzation is 



iJ^ = Ji u^v \nm^ + . . . 
and therefore the infrared divergence in the full theory amplitude is 



iA 



IR 



11 

Att' 



tanh ^ 



1 



In m,^ + . 



9' 



+ 1 

Ji u^v Inm^ + . . . 
Jiula-'vh lnm<^ + . . . , 



where Uh and Vh are two-component spinors. Note that Eq. ( |3.9|) may be written 



iA 



IR 



JiU^v [r{w) — l]lnm^ -|- . . . 



where 



and the function r{w) is given by 



w 



1 



1 + 



Vw^ - 1 



:ln 



w 



+ Vw^ - 1 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



This is the analytic continuation to the production region (negative w = v ■ v' < 0) oi 
the infrared divergence encountered in HQET (for scalar exchange) [jl2|. The first term in 
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FIG. 3. Infrared divergent one loop diagrams in NRY. The dashed hne corresponds to the 
potential scalar while the solid line is the radiation scalar. 



Eq. ( |3.9| ) is singular as /3 — > 0, corresponding to the infinite complex anomalous dimension 
found in HQET at threshold. Since it is imaginary, it does not contribute to the decay rate 
at 0{g'^). The second term in Eq. ( p.9| ) cancels in physical matrix elements with scalar 
bremsstrahlung. 

In order to be able to match onto NRY at one loop, both of these divergences must 
be reproduced in the low energy theory (the singularities higher order in (3 will only be 
reproduced when operators higher order in the velocity are included in the effective theory). 
Diagrams with both potential (fp and radiation (ppi scalars contribute to the amplitude in 
NRY (Fig. 15). The wavefunction graphs with ipp exchange vanish, while the vertex graph 
gives 

d'^k 1 



- / ^2^^, (T + A:o - MEl! + (T - A:o - Me)! + ^e) (k^ + 



— is 
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(3.13) 



where T = E — m is the fermion kinetic energy. Closing the k^ integral in the upper 
half plane, using the leading order equation of motion T = p^/2m and picking out the 
corresponding pole leaves a, d — 1 dimensional Euclidean integral, which gives 



r(2-^^) ,1 



iAp = g mJiu\a^Vh — ), , , , / dx {m {1 — x) — x p — ie 

[At:)"- ^ Jo 



= -——Jiu{(T'vh\'o^m^ + ... (3.14) 
and reproduces the first term in Eq. p.9|). The radiation scalar vertex correction is 



'-^^^ = -^^''^-"^"^"^/(^ (p-4+^^)feg ' ^'-''^ 

where the +ie^s in the fermion propagators have been dropped because these poles do not 
contribute to physical matrix elements. One can see by working with off-shell states that 
the poles from the fermion propagators only give contributions proportional to powers of 
E — 'p^ /2m which vanish by the lowest order equations of motion of NRY. Using the standard 
HQET trick of combining denominators with a dimensionful parameter, the vertex graph 
becomes 
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iAft = 8ig Jiulcr^Vh / XdX 
Jo 



1 



(27r) 



47^2 



Jiulcr'vh Inm^ + . . . . 



(3.16) 



The wave function graphs give an identical contribution, so the sum of radiation scalar 
graphs becomes 



iA 



R 



9_ 

27r- 



-Jiu\(T''vh Inm^ + . . . , 



(3.17) 



reproducing the second term in Eq. 

This illustrates that both ^pp and are required for the difference between the matrix 
elements of the external current in the full and effective theories to be infrared finite. Hav- 
ing demonstrated this, it is easier to calculate the matching conditions by regulating both 
the infrared and ultraviolet divergences in the full and effective theories with dimensional 
regularization. The matrix element in the full theory is found to be 



iA 



full 



Jiulcr'vh 1 



+ i 



3_ 

47]-2 

2 

d-A 



2 m 

+ 7£; + log 



d-A 



- ZTT + 7£; + log 



Oiv) 



(3.18) 



Regulating the theory in both the ultraviolet and infrared with dimensional regularization 
has the advantage that one-loop graphs in the nonrelativistic theory containing radiation 
scalars vanish. This is due to a cancellation of infrared and ultraviolet divergences, since 
one-loop integrals containing radiation scalars are of the form 



/ 



d'^k 
(27r)^ 



(3.19) 



which has no mass scale and so vanishes in dimensional regularization. Thus, radiation 
scalars do not contribute to the one-loop matching conditions. This does not mean, how- 
ever, that radiation scalars are irrelevant. Integrals of the form ( |3.19| ) are both ultraviolet 
and infrared divergent, with the divergent terms having the same magnitude but opposite 
signs. In the difference between the full and effective theories the infrared divergences in 
the two theories cancel, leaving an ultraviolet divergence in the effective theory. Unlike the 
infrared divergence, the ultraviolet divergence is cancelled in the low-energy theory by a 
local counterterm. 

The matrix element of the current in the effective theory with the tree-level matching is, 
in dimensional regularization. 



9' 



16np 



d-A 



irr + 'yE + log 



+ Oiv) 



(3.20) 



The difference between the two matrix elements ( p. 181 ) and ( |3.2CI| ) is analytic in the external 
momenta, as it must be to be absorbed into the coefficients of local operators in NRY. The 
matching condition for the current at one loop at a renormalization scale /i is therefore 
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2 TTl 
Ja'lpYij ( 1 - ^ +1E + log 




2 



+ 





In^^ J.V^t^^Xh, (3.21) 



where we have denoted the bare operator by the subscript 0, and the effective theory is 
renormahzed using MS. 

IV. NRQCD 

The matching of an external vector current in NRQCD, relevant for e^e~ — >■ hadrons near 
threshold, proceeds in much the same way as in the Yukawa theory of the previous section. 
Infrared divergences odd in v are reproduced in the nonrelativistic theory by potential gluon 
exchange, while infrared divergences even in v are reproduced by radiation gluon exchange. 
In this section the matching conditions for a vector current in NRQCD to order v"^ are 
calculated. The theory is regulated in both the infrared and ultraviolet in dimensional 
regularization. Cancellation in the matching conditions of terms which are not analytic in the 
external momenta is rather nontrivial at subleading order and provides a nice demonstration 
of the consistency of this approach. 

When working at subleading orders in f , there are a few subtleties which must be taken 
into account. First of all, since the three-momentum in the effective theory is 

IpI = m-i(3 = m (4.1) 
Vl 

derivatives acting on operators in the nonrelativistic theory give factors of 7/3, rather than 
f3. It is therefore more convenient to treat 7/3 as the nonrelativistic expansion parameter. 
In the rest of this paper terms of order |p|"/m" = 7"/^" will be referred to as being of order 

Secondly, Feynman diagrams in the full theory yield S matrix elements evaluated between 
relativistically normalized states, satisfying {k'\k) = 2Ek{2Tc)^5^^\k — k'). However, in 
the nonrelativistic theory defined such that the residue of the pole in the propagator is i, 
Feynman diagrams instead give S matrix elements between nonrelativistically normalized 
states, defined such that {k'\k) = 2m (2TrY S^^\k — k'). To demonstrate this, consider the 
two-point functions in the full and effective theories. Expanding the relativistic propagator 
for quarks in terms of the low-energy variables gives 

iii) + m) 2im 



p2 — TfYii {m + T)2 — p2 — m? 
2i m 

~ 2mT + T2 - p2 



^See also Ref. [|T3|, where the matching of an external electromagnetic current to nonrelativistic 
quantum mechanics (regulated in position space) was discussed. 
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where an irrelevant constant term has been dropped. This two-point function is reproduced 
in the nonrelativistic theory by a Lagrangian 



where the field operator -0^ {ip'^) annihilates (creates) a nonrelativistic particle. However, 
the residue of the pole in T — p^/2m in this theory is not i, but i{l — p^/2m + ...) = im/E. 
While this is perfectly consistent, it is preferable to remove this extra factor of m/E. This 
may be easily done, since the operator 

-"''SI'* + 3^!. 

is proportional to the equations of motion, and may therefore be removed by the field 
redefinition 

^'n^^.= (l- ^) i^'. = + Oiv% (4.5) 



However, because of this rescaling, an additional Feynman rule of y E/m for each external 
leg must be included in NRQCD, when evaluating matrix elements between relativistically 
normalized states. If this term is omitted, Feynman diagrams in NRQCD correspond to 
matrix elements between nonrelativistically normalized states. (A more careful analysis 
using the LSZ reduction formula in the nonrelativistic theory instead of rescaling arguments 
reproduces this result.) In the rest of the discussion the iph fields will be used, and matching 
conditions are calculated using nonrelativistically normalized states in the full and effective 
theories (this is the origin of the factors of yJm/E in the matching conditions presented in 
•) 

The kinetic term for the nonrelativistic fields therefore takes the usual form 

Xh 



g 



+^ (^ivV. - xi^'Xh) + 0{v'). (4.6) 

The low energy theory contains both potential {Ap = Ap^T"') and radiation (A^ = Aj^T"') 
gluons. The kinetic term for the gauge fields is, including the gauge fixing terms, 

= -\ {V^A^ - V^A^pf + I {VAlf - ^ (V • Apf 

—G^kG^^u - ^ ' + 0{g^v) , (4.7) 

where Coulomb gauge corresponds to the limit a ^ 0. The 0{gy/v) terms correspond 
to the triple-potential gluon vertex and will not be required for the one-loop matching we 



consider in this section. Also note that in Coulomb gauge there is no O(w^) correction to the 
A^p propagator. In Lorentz gauge, the gauge-fixing term for the radiation fields is instead 
~2a ("^Ai^H)^' there are additional terms bilinear in the Ap's suppressed by powers of v 
coming from the expansion of the gauge-fixing term. 

Working in the centre of mass frame pi = — p2 = p, and using the relation (easily verified 
with on-shell bispinors) 



u{pi)Yv{p2) =ul{l + 



E 



m 



ul I cr" - 



2m^ 



(4.8) 



the tree-level matching conditions for an external vector current J{x) are found to be 

J^^^M^ ^ CiOi + C2O2 + C3O3 + . . . + 0{v') 



Cl = C2 = 1 + 0(/) 

C3 = 0{g') , 



(4.9) 



where 



Oi = Jiipla-'xh 



O2 
O3 



4^2 
1 

2m2 



Ji 



Ji 



(4.10) 



and only the terms contributing to quark-antiquark production have been included. The 
operators in Eq. ( |4.1CI|) are renormalized in MS. 

The A^p coupling to heavy fermions, giving the Coulomb potential, is 0(f^^/^), 

Cc = -g {^PiA%^l,, + xlA'pXh) + 0{g'). (4.11) 

The leading corrections to this are the Darwin and spin-orbit couplings, arising at 0{v^^'^), 



D,SO 



(^T^. + xlT'^Xh) V'A 



Oa 
P 



(4.12) 



Transverse potential gluons couple through the p • A and Fermi (chromomagnetic dipole 
moment) terms at 0(f^''^), 



r 



p-A,F 



(4 (Ap • V + V ■ Ap) - xl (Ap ■ V + V ■ Ap) Xh) 
--^ {ijlcT ■ (V X Ap)^, - xio- ■ (V X Ap)xh) + Oig') 



g 

2m 

g_ 

2m 



(4.13) 



and so transverse gluon exchange and the leading relativistic corrections to Coulomb ex- 
change both contribute to potential scattering at 0{g'^v), as expected. This is also the 
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FIG. 4. One loop contributions to quark-antiquark production in QCD. 



same order as the correction to Coulomb scattering due to the correction to the fermion 
legs. This agrees with the power counting of Ref. |^ in which the matrix element of each 
of these terms is given as 0{v^), since in quarkonium v ^ g"^. Note also that the Fermi, 
Darwin, spin-orbit and relativistic kinematic corrections in the previous equations are only 
the leading pieces of the usual form of these terms 



6C 



bilinear 



+ 



i(D -gE-gE- D)z^, + xI(D ■ gE - gE ■ Ei)xh] 
+ ^ X - X iB) ■ (T^h + xI(«D X ^E - ^E X iD) ■ axh) 

+ ^ {i;i{gB ■ (T)ij, - xiigB ■ C7)xh) + 0{g^) , (4.14) 

since covariant derivatives and E both consist of two terms of differing orders in v. 

Since radiation gluons do not contribute to one-loop graphs, as was discussed in the 
previous section, their couplings are not presented here. 

Using the relations ( |4.8|) and 



2m 



(4.15) 



(still working in the centre of mass frame), the amplitude for quark-antiquark production in 
QCD from the diagrams in Fig. § may be expanded in powers of p/m: 



, f t i 



127r2 



m 

IpI 



2m^ 



mIp ■ crp'vh 1 



37r2 



+ 



+ 



3M 
2m 

P^ 



\^ + in \^E + 

vr^ + i-K {^E - 2 
(2 „ 16 



2 



d- A TTfl^ 
,2 
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FIG. 5. One loop contributions to quark-antiquark production in NRQCD. The dashed hne 
corresponds to a potential Aq gluon, the dashed gluon line to a potential Ai gluon. The shaded 
circles represent (b) the p- A vertex, (c) the Fermi vertex, (d) the Darwin vertex, (e) the relativistic 
kinematic correction to the fermion leg, and (f) O2. Implicit in both (d) and (e) are graphs with 
the same operator insertion on the antiquark line. The wavefunction graphs vanish. 



(where an overall factor oi E/m has been divided out, to convert to nonrelativistically nor- 
malized states). The amplitude has the expected 1/v singularity from Coulomb scattering, 
signalling the failure of perturbation theory close to threshold. 

The matching conditions for the current are given by the difference between Eq. ( [4.16[ ) 
and the graphs in Fig. |^ computed on-shell in NRQCD [15|. The graph in Fig. |^(a) 



corresponds to Coulomb Ap exchange. The only subtlety in evaluating this graph is that, 
to the order in which we are working, the on-shell condition in the effective theory is 

T=|--^, (4.17) 

instead of the leading order relation T = p^/2m, and this extra term must be treated 
correctly as a perturbation so as not to violate power counting. Evaluating the graph in 
Fig. |^(a) for arbitrary off-shell spinors 

'"^^ 1 + S, (4.18) 



p2 



gives 



(4.19) 

The hypergeometric function 2-^1(0, b,c]C,) may be expanded in powers of S, giving 
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(4.20) 



Evaluating this graph with the leading-order on-shell condition, 5 = 0, the 0(5'^"^) term 
vanishes as long as Iie{d) > 4. In order to have this result remain as the leading term in 
the expansion away from 5 = so that power counting is retained, the 5*^"^ term must be 
evaluated in dimensional regularization as a formal power series, f{6) = /(O) + Sf'{0) + . . .. 
In this case, each term in the expansion vanishes for sufficiently large Re((i), so the entire 
series vanishes in dimensional regularization. The final result for this graph near d = 4 is 
therefore 
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(4.21) 



which, for 6 = — p'/4m', gives 
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(4.22) 



This reproduces the 0(1/ v) term in the full amplitude. 

As discussed above, there are no graphs at 0{g'^v'^"') in NRQCD from radiation gluon 
loops. At 0{g'^v) there are contributions from the leading relativistic corrections to Coulomb 
scattering. In addition, since Coulomb exchange scales as v~^, the dressing of O2 with a 
single Ap exchange also contributes at 0{g'^v). Ap exchange contributes both via the p ■ A 
coupling (Fig. Kb)) 
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(4.23) 



and the Fermi coupling 
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Coulomb exchange is corrected by the Darwin vertex, 
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(4.25) 



while the spin-orbit coupling does not contribute. The relativistic corrections to the quark 
and antiquark propagators give 
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and finally, the one-loop correction to O2 in Fig. |^(f) gives 
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Combining these results gives 



^v4nrqcd = ciulcr'vh 
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As required, all the nonanalytic dependence on the external momentum cancels in the match- 
ing. This result is also gauge independent. 

Comparing Eqs. (|4.16| ) and ( [4.28|) gives the coefficients ci — C3 (regulating the low-energy 
theory as usual in MS) to 0{g'^): 
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C2 
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Oial 



Sin- 



(4.29) 



The result for ci reproduces the familiar short-distance correction to e^e~ — qq near thresh- 
old whereas C2 and C3 generalize this to O(f^). Note that the bare ci is finite while the 
bare C3 is divergent. This reflects the fact that there are no infrared or ultraviolet diver- 
gences in the amplitude at 0{v^) since the quarks are in a colour singlet state, and therefore 
cannot radiate a gluon at leading order in the multipole expansion. 

The major result of this section is that the nonanalytic dependence on the external 
momenta in the QCD amplitude is exactly reproduced in NRQCD. This provides a nontrivial 
check of the consistency of this approach beyond leading order. However, only for values of 
the coupling and external momenta such that ^ f <^ 1 does the tree-level matching of 
O2 and O3 dominate the two- loop matching of Oi. For scattering states closer to threshold 



(as well as for bound states) where v 



< 



a. 



ladder graphs containing potential gluons must 



be summed to all orders via the Schrodinger equation. In this case, graphs containing a 
single insertion of the tree- level matching of O2, the two- loop matching of O2 and the tree- 
level matching of Oi combined with a single higher order potential contribution are equally 
important. In this region, the one- loop matching to O2 and O3 that are presented here are 
the same order as the three loop matching to Oi. 
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V. CONCLUSIONS 



In this paper we have presented a power counting scheme for nonrelativistic effective 
theories that allows for a systematic calculation of subleading effects. A systematic v count- 
ing scheme simplifies the calculation of relativistic corrections to QCD processes such as 
quarkonium production and decay. As it is usually presented, NRQCD does not have man- 
ifest V power counting in the Lagrangian, nor is v power counting preserved by loop graphs 
either in dimensional regularization or with a momentum cutoff. While there is nothing 
in principle wrong with this, it makes calculating matching conditions somewhat awkward, 
since the matching conditions for any given operator will change by 0(1) when higher order 
operators are included in the Lagrangian. 

Velocity power counting is only preserved by loop graphs in dimensionally regulated 
NRQCD when gluons contributing to potential scattering are treated separately from on- 
shell gluons. This was accomplished in this paper by introducing two distinct gluon fields in 
the effective theory. Potential gluons propagate instantaneously and give rise to the QCD 
potential, whereas radiation gluons do not contribute to potential scattering, but correspond 
to on-shell gluons. The power counting is manifest in the Lagrangian when space, time and 
the fields are rescaled for the potential fields as discussed in Ref. [0 and for radiation fields 
as discussed in Ref. |TT|. As shown in Ref. |TT[, under this rescaling radiation fields couple 



to fermions via the multipole expansion. Separating these gluon modes realizes at the level 
of the Lagrangian the separation advocated for NRQED in Ref. 0. Under this rescaling 
V power counting is manifest in any gauge, not just Coulomb gauge, and also holds for 
non-gauge interactions. 

The infrared divergences arising in fermion-antifermion production in Yukawa theory at 
order and v'^ were shown to be reproduced in the nonrelativistic effective theory only 
when both potential and radiation scalars were included, and the matching conditions at 
that order were shown to be analytic in the external momentum. Finally, the matching 
conditions for quark-antiquark production by an external vector current were computed in 
NRQCD to 0(^2^2) 
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